Equilibrium and non-equilibrium fluctuations in a glass-forming liquid 
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Glass-forming liquids display strong fluctuations - dynamical heterogeneities - near their glass 
transition. By numerically simulating a binary Weeks- Chandler- Andersen liquid and varying both 
temperature and timescale, we investigate the probability distributions of two kinds of local fluc- 
tuations in the non-equilibrium (aging) regime and in the equilibrium regime; and find them to be 
very similar in the two regimes and across temperatures. We also observe that, when appropriately 
rescaled, the integrated dynamic susceptibility is very weakly dependent on temperature and very 
similar in both regimes. 
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The study of fluctuations in statistical mechanics has 
had a long history, starting from Einstein's seminal work 
on Brownian motion. Recent research on fluctuations 
has focused on nonequilibrium problems, such as turbu- 
lence [1], non-equilibrium steady state phenomena [2], 
single molecule experiments in biophysics [3], or the re- 
laxation of large perturbations [4] . In the present work, 
we want to use glasses as a convenient laboratory to study 
fluctuations in both equilibrium and out of equilibrium 
states; i.e. on both sides of the (kinetic) glass "transi- 
tion" . 

As experimentally observed, the glass transition is not 
a sharp thermodynamic phase transition [5], but rather 
a dynamical crossover. As the system approaches this 
crossover, the equilibration time r eq gradually increases, 
until it eventually becomes longer than the representa- 
tive laboratory timescale t\ a t, [6, 7]. Beyond the glass 
transition, the system cannot equilibrate in laboratory 
timcscales. This leads to the emergence of new phe- 
nomena, such as physical aging, i.e. the breakdown of 
time translation invariance (TTI). Besides the slowdown 
of the dynamics, also dynamical heterogeneities [7, 8], i.e. 
strong local fluctuations in the properties of glass-forming 
liquids, are observed both in the equilibrium regime [9, 
10, 11] and in the aging regime [12, 13, 14, 15, 16, 17]. 
It is not clear, though, to what degree the properties of 
those fluctuations depend on whether they are observed 
in one or the other regime. In the present work, we ad- 
dress this question in the context of a detailed molecular 
dynamics simulation of a simple glass model probing the 
aging regime, the equilibrium regime, and the crossover 
regime between the two. 

A theoretical framework based on the presence of a 
Goldstone mode in the aging dynamics [18] that gives rise 
to local fluctuations in the age of the sample [19, 23] pre- 
dicts that probability distributions of local fluctuations 
in the aging regime are approximately independent of the 
waiting time t w at fixed values of the two-time global cor- 
relation C g iobai(i, t w ). This prediction, initially proposed 
for spin glasses [19], has been found to apply also to off- 



lattice models of aging structural glasses [15, 17] (where 
Cgiobai is the self part of the intermediate scattering func- 
tion: C g iobai(M™) = ^E^Liexp^q-faW - r j(*«>))) )■ 
Additionally, it has been found that the dynamic spatial 
correlations in an aging structural glass display a simple 
scaling behavior as a function of C g i bai [16, 17]. In what 
follows we address the question of whether these same 
scaling behaviors extend to the equilibrium regime of an 
off-lattice structural glass model. 

In order to observe both equilibrium and aging phe- 
nomena, one could either reduce the temperature T of the 
thermal reservoir that the system is coupled to beyond 
the glass transition temperature T g at a fixed laboratory 
timescale ti a b, or reduce ii a b until it becomes shorter than 
T eq , at fixed reservoir temperature T. Here we use both 
approaches. To observe genuinely glassy dynamics, we 
make sure that T is low enough that two-step relaxation 
is observed, and that timescales are much longer than the 
timescale of the first step. 

We simulate a 80 : 20 binary mixture of par- 
ticles of mass Ma = Mb — M, in three di- 
mensions, at a number density p = 1.204, in- 
teracting via short-range, purely repulsive Weeks- 
Chandler- Andersen (WCA) potentials [20]: V a .p = 



^[(^Y-(^y + l\9(r^ oa > - r), where 

a,/3 G {A, B} denotes the particle type, e Qj( 3 and a a jj 
are taken as in [21, 22], 9{x) is the Heaviside theta 

function, and r = 2 1 ^ 6 a a ,p « 1.1225 makes 

V(r) and dV/dr continuous for all r > 0. We take 
ks = 1.0, and the units of length, energy and time as 

1/2 

<JAA, tAA and t lj = (a A A 2 M/48e A A) respectively, 
where tlj represents roughly the timescale for vibra- 
tional motion in the system. Initially the system is equi- 
librated at a high temperature: Tj = 5.0, and then in- 
stantaneously quenched to a final temperature Tf < T;. 
The top panel of Fig. 1 illustrates the temperature and 
time regimes for which we show results. We take three 
different T/'s. For T f = 0.236 « 0.9Tmct, only the ag- 
ing regime is observed within our simulation time frame. 
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FIG. 1: Top panel: Timescales as functions of the temper- 
ature. "+" symbols: (temperature, waiting time) pairs for 
which we show our results. Triangles: a relaxation times for 
T = 0.29, ■ ■ ■ , 0.40. Full line: fit of r a to r Q = a(T-T MCT )" 7 . 
Dotted line: r eq , estimated as r eq ~ 10r a (see text). Vertical 
dash-dotted line: Imct- Bottom panel: Probability distri- 
bution p(C r ) at T = 0.29, when C global (t, t w ) = 0.1,0.3,0.5, 
for t w /r a — 0.31,0.79 (aging) and t w /r a = 11.5,15.8 (equi- 
librium). The two equilibrium plots overlap perfectly. The 
coarse graining region contains on average 6.6 particles. 



FIG. 2: Rescaled probability distributions o~cp(C T ) as func- 
tions of the normalized fluctuation (C r — Global )/°"c in the 
two-time correlation, shown at t w /r a — 0.79 (aging regime) 
and t w /r a — 11.5 (equilibrium regime), for temperatures 
T = 0.4 (triangles), T = 0.29 (circles) and T = 0.236 
(squares), with Cgi^^t, t w ) = 0.3. Generalized Gumbel fits 
to the distributions for T = 0.4 are shown with thin lines. Top 
panel: Coarse graining region containing on average 6.6 par- 
ticles/box. Bottom panel: Coarse graining region containing 
on average 30.5 particles/box. 



For T f = 0.29 w 1.1T MC t and T f = 0.4 « 1.5T MC t, 
we observe both the aging and equilibrium regimes. 
Here, T MC t denotes the Mode Coupling Theory (MCT) 
critical temperature. We define r a by the condition 
limt^-xxj C(t w +T a ,t w ) = 1/e and find it for the tem- 
peratures T = 0.29, 0.30- ■ -0.4. We observe that for 
2.9 < T < 3.4, the equilibrium C(t,t w ) displays two 
step relaxation and also scaling behavior as predicted by 
MCT. Following [22], we fit r a = (T - T MC t)" 7 in that 
range, and find T MC t = 0.263 ± 0.01 and 7 = 2.1 ± 0.7. 

The time of the quench is taken as the origin of times 
t = 0. In our simulations, t\ a b is replaced by t w . We find 
that aging effects are strongest for t w /r a < 2. For 2 < 
tw/ T a ^ 10 the aging effects gradually become weaker 
with increasing t w until they disappear for t w /r a ~ 10. 
We thus estimate r oq ~ 10t q . We show results for a 1, 000 
particle system, with 5000, 9000 and 10, 000 independent 
thermal histories when T f = 0.236, T f = 0.29, and Tj = 
0.4 respectively. 

We present results for the probability distribu- 



tions of obscrvablcs which probe local fluctuations 
in small regions of the system: the local coarse 
grained two-time correlation function [15] C r (t,t w ) = 

jvfk) E rj (t ro )es r cos (q ■ ( r iW - r i (**>)))' and thc P^- 
ticle displacements Axj(t, t w ) = Xj(t) — Xj(t w ) along one 
direction [10, 15, 24]. Here B r represents a small coarse 
graining box around thc point r in the system and N(B r ) 
is the number of particles present at the waiting time t w 
in the box B r . We choose a value of q that corresponds 
to the main peak in the structure factor S(q) of the sys- 
tem, q = 7.2. In order to probe the spatial correlations of 
the fluctuations, we also consider the generalized 4-point 
density susceptibility X4 = J d 3 r g4(r,t,t w ) [11, 16, 25], 
where gi(r,t,t w ) is a 4-point (2-timc, 2-position) corre- 
lation function [11, 16]. 

In the bottom panel of Fig. 1 we show the probability 
distributions p(C T ) for C g i b a l (t,t w ) = 0.1,0.3,0.5. We 
observe that, for fixed value of Cgi Y> a \(t,t w ), the proba- 
bility distributions are approximately invariant between 
the two regimes. These results extend those found for the 
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FIG. 3: p(Ax) for C g i b a i(i, tw) = 0.5, plotted with a 
logarithmic vertical axis to emphasize the tails of the dis- 
tributions. Fits to the tails by a nonlinear exponential 
form are shown with full lines. Top panel: T — 0.29. 
t w — 0.31 r a , 0.79 t q are in the aging regime; whilst t w = 
11.5 r a , 15.8 T a , which give perfectly overlapping results, are 
in the equilibrium regime. Bottom Panel: T = 0.23. Only 
the aging regime is accessible; t m = 10 2 , 10 3 , 10 4 . 



aging regime of a binary Lennard- Jones (LJ) glass [15]. 

Fig. 2 shows the rescaled probability distributions 
<Jcp{C r ) versus the normalized fluctuation (C r — 
CgiobaO/cc m the one-point two-time correlator, for 
coarse graining regions of two sizes, for the temperatures 
T = 0.23, T = 0.29 and T = 0.4. We observe a very 
good collapse of the data for different t w and different 
temperatures for the smaller coarse graining size, and a 
slightly less good collapse for the larger coarse graining 
size. This suggests that the weak dependence on t w may 
be due to the time dependence of the dynamic corre- 
lation length [15]. In [23] it has been argued that the 
probability distribution of the normalized fluctuations 
(C r — CgiobaO/cc can be described well by a generalized 
Gumbcl distribution $Gumbei [!]■ Fig- 2 shows that fits 
of our data to ^Gumbcl become better for a larger coarse 
graining region, i.e. when averaging affects have started 
to modify the distribution [1]. 

In Fig. 3 we show the probability distribution p(Ax) 
of the particle displacements, with C g i bai(£, t w ) = 0.5, 
for both Tf = 0.29 > T MCT (top panel) and T f = 0.23 < 
Tmct (bottom panel) . We use semilog plots to emphasize 



FIG. 4: Top panel: X4.{t,t w ) as a function of t — t w for 
T — 0.29. Bottom panel: \i divided by its maximum 
value X4( max )i plotted as a function of 1 — C, with C = 
Cgiobal(Mi") for T = 0.4 (triangles), T = 0.29 (circles) and 
T = 0.23 (squares). 



the tails of the distributions. The distributions are very 
similar to those determined by confocal microscopy for 
an equilibrium stcrically stabilized colloidal liquid near 
its glass transition [10], and also to the ones found in sim- 
ulations of an aging binary LJ glass [15]; but markedly 
different from those found by confocal microscopy in an 
attracting colloid [24]. For small to moderate values of 
Ax, the data for different t w and temperatures collapse 
very well with each other and with a common gaussian 
fit, as long as C g i bai is constant (not shown). We find, 
however, that for fixed C g i bai(i, t w ), the tails of the dis- 
tribution become wider for increasing t w . The tails can 
be fit in the region |Ax| > 0.5 by a nonlinear exponential 
form p(Ax) w N exp(— \A.x/a\^) 1 with j3 ~ 1 decreas- 
ing monotonically with increased t w [15]. For the case of 
T > Imct both the tails and the value of f3 gradually 
saturate as the system approaches equilibrium. 

To quantify the spatial correlation of dynamical fluc- 
tuations, the 4-point density susceptibility Xi{t)tw) [11, 
16, 25] and similar quantities have been used both 
in numerical simulations [11, 16, 25] and in experi- 
ments [26, 27, 28, 29]. The top panel of Fig. 4 shows 
Xi(t,t w ) as a function of t — t w at Tt = 0.29, for vari- 
ous t w . We have found that X&{t->t w ) has a peak whose 
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height and location grow with t w as long as the system 
is aging. This behavior as t w is increased is observed ex- 
perimentally in a coarsening foam [26] , and analogous to 
the behavior of supercooled liquids as the temperature is 
reduced [27] and of granular systems as the area fraction 
is increased [28] . In the bottom panel of Fig. 4 we test 
for a possible scaling behavior with C = C g i bai , by plot- 
ting the rcscalcd quantity X4/X4( max ) a s a function of 
1 - C, for Tf = 0.23, 0.29, 0.4. An exact collapse would 
implicate the factorization Xiit, t w ) = xt{t w )4>(C(t, tw)), 
where xl{t w ) = X<i(t,tw)\c(t,t w )=i/e is a rescaling fac- 
tor that depends only on t w until it saturates when the 
system equilibriates; and 4>(C(t,t w )) is a scaling func- 
tion which depends on times only through the value of 
the intermediate scattering function C(t, t w ). We observe 
that the curves approximately collapse into two groups, 
one for T > Tmct and another for T < Tmct, with a 
slight horizontal shift between the two groups. Curves for 
the same temperature approximately collapse, although 
a slight systematic variation with t w remains. 

In summary, we have analyzed the probability dis- 
tributions of local two-time observables and the spatial 
correlations of fluctuations in the equilibrium and aging 
regimes in a simple glass model. We do find some dif- 
ferences between the two regimes. For example, the tails 
of the probability distributions of one-particle displace- 
ments p(Ax) become slightly wider as t w grows at con- 
stant CgiobaitA t w ). Also, the rescaled generalized sus- 
ceptibility X4/X4( max )j when plotted as a function of 
1 — Cgiobai, shows slight systematic dependences on t w 
and on T. However, those differences are minor, and 



the similarities arc quite striking. In general terms, the 
simple scaling behaviors previously observed in the ag- 
ing regime [15, 16] extend to the equilibrium regime. 
On one hand, both the probability distributions of local 
two-time correlations and the probability distributions 
of one-particle displacements are approximately invari- 
ant between the two regimes as long as the total cor- 
relation C g iobai(i, t w ) is kept constant. On the other 
hand, the rescaled dynamic susceptibility X4/x4(max), 
when plotted against 1 — C g i bai(A t w ), is also approxi- 
mately invariant between the two regimes. Additionally, 
we find that the probability distributions p(C r ) of local 
two-time correlations can be fitted by generalized Gum- 
bel forms [1, 23], as long as the coarse graining size is 
relatively large. Finally, the scaling functions describing 
the distribution of local two-time correlations and the dy- 
namic susceptibilities are very weakly dependent on tem- 
perature across the range from T < Tmct to T > Tmct, 
and qualitatively similar to those found in an aging LJ 
glass [15, 16]. 
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